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Abstract 

Techniques  for  implementation  of  a Nyquist  stability  result  for  a linear  time 
invariant  multivariable  feedback  system  are  described.  The  approach  is  based  on 
continuation  methods  for  computing  the  system's  eigenvalue  loci. 


I.  INTRODUCTION 

The  classical  Nyquist  stability  criterion  for 
single-input  single-output,  linear  time-invariant 
feedback  systems  has  only  recently  been  generalized 
to  multivariable  feedback  systems  [1,2].  Stability 
theorems  are  expressed  in  terms  of  the  eigenvalue 
loci  of  the  open  loop  transfer  function  G(s)  of  the 
system.  In  particular  if  G ( s ) is  stable,  i.e., 

G(s)  has  no  poles  in  the  right  halF  of  the  s-plane 
or  on  the  ju-axis,  then  a linear  time-invariant 
multivariable  feedback  system  with  n inputs  and  n 
outputs  is  stable  if  and  only  if  its  generalized 
Nyquist  plots  (union  of  eigenvalue  loci)  does  not 
pass  through  or  encircle  the  (-1,  0)  point  [1]. 

In  order  to  apply  the  multivariable  Nyquist  criter- 
ion, it  is  thus  necessary  to  compute  the  eigen- 
value loci  as  a function  of  frequency.  For  a given 
frequency,  the  eigenvalues  can  be  calculated  by 
using  classical  techniques.  Since  the  eigenvalues 
are  functions  of  frequency,  normally  one  would  have 
to  repeat  the  entire  computational  procedure  for 
each  frequency.  In  the  actual  stability  analysis, 
this  repetition  is  however,  impractical.  Our 
approach  to  the  stability  analysis  of  multivariable 
feedback  systems  is  based  on  continuation  methods. 


The  basic  idea  of  all  continuation  methods  is  to 
convert  the  solution  of  a parameterized  family  of 
algebraic  problems  into  the  solution  of  a differ- 
ential equation.  Then  if  one  can  find  the  solution 
of  an  initial  problem  by  using  classical  methods  the 
solutions  to  the  other  problems  can  be  obtained  by 
integrating  the  associated  differential  equation 
with  the  initial  solution  as  an  initial  condition. 

II.  EIGENVECTOR  APPROACH 

Our  first  method  is  based  on  the  approach  described 
by  Faddeev  and  Fadeeva  [3]  and  Van  Ness  et.  al.  [4]. 
A differential  equation  is  written  with  the  eigen- 
values as  dependent  variables  and  the  frequency  as 
variable  parameter.  We  then  compute  a set  of 
initial  eigenvalues  by  classical  analysis  techniques 
and  integrate  the  resulted  differential  equation  to 
obtain  the  required  eigenvalues  for  each  frequency. 
The  eigenvalues  of  G(jw)  and  their  complex 
conjugates  X.(u)  satisfy 

G(jj)X.(u)  = X1(o.)Xi(u.)  1=1.2 n (1) 

and 

G*(j<o)Vi(w)  = A.(u.)Vi(co)  1=1,2 n (2) 

where  X.(m)  and  V.(u>)  are  the  corresponding  eigen- 

* ir 

vectors  of  X.(u)  and  X.(w)  respectively,  and  G (jw) 
is  the  complex  conjugate  transpose  matrix  of  G(j-). 
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We  differentiate  (1)  with  respect  to  u to  yield 


dx.  <5^  Xs,  V.> 

7T%;m  uV  . 1 = 1.2 n.  (3) 


= -f(t)  f (0)  = f[x(u0)]  = 0 


du  <Xit  V.> 


The  differential  equations  involving  X.  and 


u(0)  = uQ. 


are  obtained  as 


~3f 

af 

n 

ax 

3a> 

- £ aiiXi*  *=  ^ * 2»  •••*  n* 

j=l  J J 

(4) 

J= 

3u 

da) 

ax 

du> 

Assuming  the  nonsingularity  of  the  Jacobian  Matrix 

[“af  af-]  Ilf  af"l 

ax  3u  ax 


XT  = j=l 


eijV  1 = 2’ 


xi*  V 


p. . = 0,  p.  . = -rr — v — 
11  1J  <Vj,  Xj> 


i 1 j. 


Starting  with  a set  of  predetermined  initial  con- 
ditions X^O)  = X.Q,  X1(0)  = X.0  and  V.(0)  = Vio 
for  i = 1 , 2,  . . . , n,  we  integrate  (3),  (4)  and 
(5)  to  obtain  the  required  eignevalues  for  each 
frequency.  The  eigenvalue  loci  are  computed  in  a 
continuous  manner  by  numerical  integration. 

III.  JAC03IAN  METHOD 

For  an  nth  order  system,  the  above  algorithm  re- 
quires the  numerical  integration  of  a set  of  3n 
equations  and  the  computation  of  two  sets  of  un- 
wanted variables— namely  the  eigenvectors  X^  and 


\ly  These  disadvantage,  can  easily  be  avoided  if 


in  the  x-u>  space  the  algorithm  (9)  reduces  to 

rdxi  r,i  r„  i 


K = j-1 


It  is  seen  from  the  solution  of  (9) 


f(t)  = 0e_t  = 0 


cii  = + t. 


that  for  any  admissible  pair  of  uo  and  x(uQ)  satis- 


the  characteristic  equation  for  the  multivariable 
feedback  system  can  be  predetermined.  A much  sim- 
pler method  can  be  formulated  based  on  the  approach 
for  finding  multiple  solutions  for  a nonlinear 
equation  developed  by  Chao  et.  al.  [5]. 

Let  the  characteristic  equation  of  G(jw)  be  given 
by  an  nth  order  polynomial  in  eigenvalue  ).  with 
complex  coefficients 

f[x(u)]  = | XI  - G(jw) | = 0.  (8) 

Instead  of  solving  (8)  directly  for  each  frequency, 
we  consider  two  simultaneous  differential  equa- 
tions of  the  form 


fying  (8),  the  corresponding  trajectory  will  remain 
on  the  solution  curve  f-0  as  u changes.  The  + or  - 
sign  is  chosen  depending  on  whether  one  would  like 
to  increase  or  decrease  t>.  Equation  (11)  may  now 
be  solved  by  any  numerical  integration  techniques 
and  the  eigenvalue  loci  can  be  traced  automatically 
by  integrating  only  a second  order  differential 
system. 

IV.  EXAMPLE 

To  illustrate  the  approaches  presented,  consider  a 
linear  time-invariant,  multivariable  feedback 
system  with  open  loop  transfer  function  character- 
ized by 


for  which  the  characteristic  equation  is  given  by 

o l/.l.  , ic  2 i. 


f[x(u)]  = X2-8X  + -IC~-  ^ =0.  (14) 

1 + u 1 + u 


The  generalized  Nyquist  plots  shown  in  the  accom- 
panied figure  for  the  cases  where  k-16  and  36  are 
obtained  by  applying  the  eigenvector  approach  where 


■ 
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as  in  the  critical  case,  k=25,  the  Jacobian  method 
has  been  used. 

In  all  three  cases,  the  equations  are  integrated 
using  Euler's  method  with  a step  size  of  0.01.  It 
is  seen  from  the  figure  that  the  system  is  stable 
for  k<25  since  the  generalized  Hyquist  plots  do 
not  encircle  -1  point. 


Im  X 
T2.4 
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